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Summary. 

This  report  presents  experimental  data,  mode  of  operation  and  construction  details 
of  advanced  plasma  focus  machines  (APF)  for  different  optimal  levels  of  discharge-feeding 
energy  W,  in  particular  for  APF-200  (W  <  200  kJ)  and  APF-50  (W  <  50  kj).  The 
function  of  these  APF  systems  was  to  determine,  along  with  the  data  of  smaller  machines, 
the  scaling  laws  of  the  emission  (fluence)  of  ion  and  ion-cluster  beams  as  a  function  of  W, 
ejected  from  the  self-field-dominated  plasma  of  the  APF  pinch.  Typical  ion  spectra  from  a 
Thomson  (parabola)  spectrometer  in  the  80°  direction  from  the  electrode/pinch  axis  are 
also  included  here. 
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1.  Introduction. 

The  primary  objective  of  the  research  program  was  the  determination  of  the 
scaling  laws  of  the  emission  of  ion  clusters,  ion  superclusters  and  of  self-field-dominated 
fragments  of  the  current  sheet  ejected  from  plasma  focus  (PF)  discharges,  as  a  function  of 
the  discharge-powering  energy  W.  The  program  required  the  construction  and  operation 
of  several  PF  machines  of  different  magnitude.  The  tasks  completed  during  the  first  year 
of  the  program  are; 

1].  Monitoring  and  evaluation,  via  statistical  analysis  of  experimental  data,  of  the 
performance  of  the  Advanced  Plasma  Focus  system  APF-50  (W  <  50  kJ).  This  system  is 
characterized  by 

(a)  Choice  of  power  level,  construction  parameters  and  filling  pressure  such  as  to 
achieve  optimization  for  operating  within  the  interval  15  kJ  <  W  <  50  kJ  of  the  energy  W 
of  the  powering  capacitor  bank  charged  at  a  voltage  Vq,  23  kV  <  Vq  <  40  kV. 

(b)  Two-module  capacitor  bank  with  five  capacitors  (Configuration  1,  inductance 
of  each  module  L„j,  =  6  nH),  or  four  capacitors  (Configuration  2  ,  inductance  of  each 
module  Ln,^,  =  7.5  nH,  energy  of  two  modules  12  kJ<  W  <  40  kJ,  at  23  kV<  Vq  <  42  kV) 
in  each  module;  each  capacitor  has  5.6  pF. 

(c)  Low  jitter  (<  10  ns)  switch-triggering  system,  with  two  rail-gap  switches. 
Maxwell  model  No.  40200  each  of  inductance  L^^  =  20  nH.  Construction  details  of  this 
APF  machine,  with  radius  of  the  two  coaxial  electrodes,  respectively  R  =  64  mm  (inner 
radius  of  outer  electrode,  cathode  )  and  r  =  26  mm  (outer  radius  of  center  electrode, 
anode),  such  as  electrode  length,  inductance  of  the  complete  electric  APF  circuit  including 
the  powering  capacitor  bank  and  basic  diagnostics  probes,  are  reported  in  the  Appendix  1 

Optimization  means  here,  as  usual,  that  the  energy  densification  process  in  the 
pinch  provides  a  maximum  of  energy  density  in  the  pinched  plasma,  as  compared  with 
what  can  be  achieved  with  an  alternative  choice  of  the  parameters  listed  in  (a),  in  the  same 
energy  interval  W.  We  used  the  neutron  yield  per  discharge,  Y,  from 
D  +  D  nuclear  reactions,  with  deuterium  filling  of  the  APF  discharge  chamber,  as  a 
convenient  figure  of  merit  to  assess  the  energy  density  level  in  the  pinch  .  Simultaneously 
to  the  D  +  D  reactions,  a  more  general  method  for  determining  the  energy  distribution  of 
the  accelerated  ions  (trapped  in,  and  ejected  from,  the  pinch)  was  set-up  and  perfected  for 
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application  to  all  APF  systems  of  the  Laboratory.  This  general  method  is  based  on 
nuclear  activation  of  solid  targets  and  of  gaseous  targets  (implemented  with  doping  of  the 
filling  gas  of  the  APF  discharge  chamber  with  high-Z-nuclei )  and  on  a  variety  of  high-Z 
nuclear  reactions,  with  different  values  of  the  energy  threshold  of  the  reaction,  as  reported 
in  the  publication  [1], 

2] ,  Completion  of  the  construction  and  operation  of  the  Advanced  Plasma  Focus 

system  APF-200  (W  <  200  kJ)  initiated  under  a  previous  AFOSR  Grant  (No.  AFOSR  - 
91-0158),  characterized  by 

(a)  APF  optimization  via  a  choice  of  the  construction  parameters  such  as  to 
permit  operation  at  different  powering  energy  levels  W  within  a  wide  interval 

70  kJ  <  W  <  200  kJ,  with  charging  voltage  Vq,  25  kV  <  Vg  <  42  kV  of  a  five-module 
capacitor  bank  (Configuration  1).  Four  modules  were  assembled  (Configuration  2)  with 
twenty  capacitors  in  each  module  of  inductance  =  1.5  nH  (2.8  pF,  30  nH  capacitors), 
with  the  readily-available  option  of  connecting  to  the  APF  fitting  header  a  fifth  module 
with  additional  20  capacitors  on  the  back  side  of  APF-200,  depending  on  experimental 
requirements.  This  adaptability  of  the  APF  design  to  satisfy  a  variety  of  experimental 
requirements  is  one  of  the  new  features  of  APF,  including  field  emitters  for  controlling 
breakdown  conditions  on  the  pyrex  insulator  sleeve  [specifically  a  circular  knife  edge 
(KE)  at  the  breech  side  of  the  electrodes],  which  differentiate  APF  from  PF  machines. 
Each  module  was  singularly  tested  with  a  rail  gap  switch  [as  for  APF-50]  and  a  pulse 
trigger  generator  [Maxwell  mod.  No.  40295]  capable  of  simultaneously  triggering  up  to 
ten  Maxwell  rail  gap  switches. 

(b)  Pentagon-shaped  APF  header  (formed  of  power-transmission  back  plates, 
electrodes  and  discharge  chamber)  connected  to  four  modules,  each  of  inductance  1.5  nH. 
Ultimately,  the  operation  of  APF-200  was  limited  to  the  simultaneous  discharge  of  only 
two  modules  out  of  five,  because  only  two  Maxwell  switches  instead  of  the  planned  five 
switches,  were  available  untill  the  end  of  the  first  year  of  the  program.  The  limitation  to 
two  switches  was  determined  by  circumstances  unforeseeable  at  an  early  stage  of  the 
research  program  including,  but  not  limited  to,  the  termination  of  the  laboratory  activity 
after  the  first  year,  instead  of  the  programmed  three-year  activity. 

(c)  Same  low-jitter  switches  and  triggering  system  as  for  APF-50. 

The  data  obtained  under  the  two-module  limitation  on  W  (35  kJ  <  W  <  100  kJ  at 
25  kV  <  Vg  <  42  kV)  are  sufficient  for  a  comparison  of  the  good  performance  of  APF- 
200  with  that  of  APF-50  at  the  same  level  of  peak  amperian  current  I  on  the  electrodes 
(I  ~  1  MA),  respectively  of  radius  R  =  97  mm  (outer  electrode,  cathode)  and  r  =  55  mm 
(inner  electrode,  anode). 

3] .  Same  diagnostics  and  set-up  for  statistical  analysis  of  data  as  for  APF-50. 

The  data  obtained  from  APF-50  have  confirmed  the  existence  of  a  particular 
(oscillatory)  mode  of  propagation  of  the  plasma  current  sheet  between  the  electrodes, 
until  the  last  stage  of  the  current  sheet  implosion  on  the  electrode  axis.  Characteristics  of 
this  oscillatory  mode  of  propagation  are  reported  in  the  publication  [2]. 
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Contributions  to  the  research  program  and  to  this  report  were  provided  also  by 
Prof.  C.  M.  Luo,  Visiting  Scientists  from  Tsinghua  University,  Beijing,  Dr.  K.  Melzacki, 
Visiting  Scientist  from  KM  Physics,  Troy,  NY,  Mr.  D.  Goldstein  and  Ms.  L.  Schneider, 
Graduate  Students  with  the  support  of  a  three-year  AFOSR-AASERT  Grant,  and  Mr. 
J.R.  Brzosko  under  a  Senior  Student  Project  in  Physics.  Mr.  M.  Kania  provided  essential 
technical  assistance  in  all  activities  of  the  Laboratory. 


2.  Performance  of  APF-200. 

Flexible  copper  plates  of  length  suitable  for  matching  the  relative  jitter  of  the 
switches,  connect  switches  and  header.  If  compared  to  coaxial  cable  for  power 
transmission,  these  flexible  plates  permit  the  use  of  higher  charging- voltage  levels  Vg  in 
air  (no  insulating  oil  was  used)  without  flash-over  at  switch  and  header  connections. 

A  Rogowski  coil  for  monitoring  the  electrode  current  I(t)  as  a  function  of  time 
was  located  between  the  two  coaxial-electrode  back  plates  of  the  header.  A  voltage- 
divider  probe  providing  the  voltage  signal  V(t),  from  the  breech  side  of  the  electrodes, 
was  connected  to  the  power  transmission  plate  attached  to  the  base  of  the  hot  central 
electrode  (anode). 

The  corresponding  two  signals  I(t),  V(t)  were  used  to  determine  the  period  T  of 
the  complete  electric  circuit  of  APF-200  [including  the  powering  capacitor  bank]  and  the 
quality  of  the  pinch  from  the  dips  of  the  I  signals  and  the  voltage  signal  spikes.  These 
voltage  spikes  occur  during  the  surge  of  anomalous  (i.e.,  non-collisional)  resistivity  in  the 
plasma  pinch  region,  during  the  pinch  disintegration  stage  and  may  have  an  amplitude 
greater  by  a  factor  >  5  than  the  maximum  charging  voltage  Vg  of  the  capacitor  bank.  The 
limitation  to  Vg  <  42  kV  avoids  flash-over  between  power  transmission  plate,  at  the 
capacitor  bushing  and  at  the  connection  between  power-transmission  plate  with  header 
and  with  switches. 

The  neutron  yield  per  discharge  Y,  from  D  +  D  fusion  reactions,  was  also  used  for 
monitoring  the  quality  of  the  energy  densification  process  in  the  pinch. 

The  general  view  of  APF-200  is  shown  in  Fig.  1 .  The  typical  dl(t)/dt  and 
I(t)  ~  /g‘  V(t')  dt'  signal  are  shown  in  Fig.  2.  The  signals  are  carried  from  the  APF  header 

to  the  monitoring  system  inside  a  noise-shielding  Faraday  cage,  via  50  Q  coaxial  cables 
with  solid-copper  double  shielding.  The  observed  period  for  the  two-modules  operation 
was  T  =  1 5.5  ±  0. 1  ps.  The  corresponding  effective  value  of  the  inductance  L  of  the 
complete  circuit,  including  the  plasma-current  sheet  propagating  between  the  electrodes, 
was  obtained  from  L  ~  (T/27i)^C  '  =  55  ±  1  nH,  where  C  is,  essentially,  the  capacitance  Cg 
of  the  utilized  two  modules,  i.e.,  C  ~  Cg  =  112  pF.  The  same  considerations  on  L,  on 
current  sheet  structure  and  multiplicity,  and  on  diagnostic  methods  presented  for  APF-50 
in  Appendix  1,  apply  to  APF-200. 

The  neutron  yield  Y  of  the  two-module  system  is  reported  in  the  histogram  of  Fig.  3. 

The  operation  conditions  of  APF-200  reported  here,  Vg  =  25  kV,  W  =  35  kJ, 
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Dj-filling  pressure  of  the  discharge  chamber  p  =  5.5  Torr,  electrode  length  to  have  the 
pinch  maximum-compression  phase  at  a  time  of  peak  electrode  current  [t  =  T/4],  are  such 
that  the  corresponding  peak  electrode  current  I  is  the  same  as  for  APF-50  with  Vq  =  28 

kV,  W  =  1 7.5  kJ,  Co  =  45  pF  (Configuration  2),  Vq  =  25  kV,  W  =  17.5  kJ,  Cq  =  56  pF, 
(Configuration  1),  specifically,  I  ~  Vo/(L/Co)'^^  ~  1  MA.  In  these  estimates,  unless 
otherwise  stated,  we  always  use  the  effective  value  L  of  the  total  inductance  of  the  circuit 
obtained  from  the  measured  value  of  T,  i.e.,  from  the  dynamic  conditions  of  the  complete 
circuit  [including  the  propagating  current  sheet],  rather  than  the  static  value  Lq  of  the 
circuit  (in  APF-50,  Lg  =  40  nH  for  Configuration  1,  Lg  =  41  nH  for  Configuration  2). 

By  definition,  Lg  is  measured  from  the  period  of  the  complete  circuit  after  shorting  the 
electrodes  at  their  muzzle  side,  without  the  plasma  current  sheet  between  the  electrodes. 

The  primary  results  are: 

(1)  The  observed  maximum  value,  Max(Y),  of  the  yield  Y  of  APF-200,  was 
Max(Y)  =  6  X  10^  of  ~  2.5  MeV  neutrons  under  the  limiting  current  condition  indicated 
above,  from  a  series  of  261  discharges  under  essentially  an  identical  choice  also  for  all 
other  operation  conditions.  This  yield  is  equal  to  the  average  value  Y^^  of  Y  for  fully- 
optimized  operation  conditions  of  APF-50  with  Configuration  2. 

(2)  No  appreciable  variation  AY^^,  of  Y^^  was  observed  by  increasing  W  from 
35  kJ  to  50  kJ  (Vg  =  30  kV),  with  a  corresponding  increase  of  I  from  1  MA  to  1.35  MA, 

i.e.. 


AY,yAW~AY,yAI~0, 


with  Y,^~Max(Y)/14,  AWAV~43%,  AI/I~35%.  This  indicates  that  APF-200  under 
the  two-module  restriction  on  the  current  and  on  Lg  operates  in  a  highly-unstable  regime 
of  current  distribution  in  the  interelectrode  gap  and,  consistently,  of  the  energy 
densification  process  in  the  pinch. 

Leading  causes  of  this  unstable  regime  is  the  relatively-low  current  density  on  the 
electrodes  [  ~  I(t)  /27i:r]  and  of  the  relatively-high  value  of  the  inductance  of  the  two- 
module  system  (module  capacitors,  switches,  power-transmission  plates), 

+  (L,„j|  +  Lj,^^)/2,  where  L,,jr  —  Lmdi  Uw 's  the  inductance  of  the  two  back  plates  of 
the  header,  [not  well  matching  the  the  coaxial-electrode  load]  as  compared  to  the 
inductance  of  the  optimizing  four-module  system. 

This  high  inductance  in  the  up-stream  part  of  the  electric  circuit  is  the  cause  of 
multiple  restrikes  on  the  APF-breech  insulator  with  a  broadening  of  the  current  density 
distribution  in  the  plasma.  This  point  is  confirmed  by  the  systematic  comparison  of  the 
two-module  APF-200  data  with  the  data  from  the  optimized  two-module 
APF-50,  Configurations  1  and  2,  tested  in  parallel  experiments. 
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1.  Performance  of  APF-50. 


Typical  diagnostic  procedures  available  during  the  operation  of  all  our  APF  and 
PF  machines  are  reported  in  Appendix  1,  where  magnetic  probe  data  and  neutron  yield 
are  reported  for  Configuration  1 . 

Fig.  4  reports  the  general  view  of  APF-50,  Configuration  1  (two  modules  with 
five  5.6  |iF  capacitors  each;  fixed  facility;  stainless  steel  discharge  chambers  with  eight 
ports)  and  Configuration  2  (two  modules  with  four  5.6  pF  capacitors  each;  mobile 
facility,  mounted  on  wheeled  platform). 

The  neutron  yield  Y  of  300  discharges  under  identical  conditions,  with  typical 
fluctuations  of  Y  from  discharge  to  discharge  is  reported  in  Fig.  5.  A  similar  distribution 
was  obtained  for  Configuration  1 . 

The  discharge  distribution  of  the  optimized  APF-50  has  a  relatively-narrow  peak 
on  a  0. 1  <  YA"av  <3.0  interval,  with  Ygv( APF-50)  ~  Max[Y(non-optim.APF-200)]  = 

6  X  109  neutrons  per  pulse  and  Max  [Y(APF-50)]  =  1.7  x  10>o  neutrons  per  pulse  from 
D  +  D  reactions  ( ~  2.5  MeV  neutrons).  A  silver  activation  counter  (Los  Alamos  type 
with  four  Geigers)  monitors  the  neutron  yield  per  pulse.  Standard  scintillations  detectors 
(Pilot-U  and  NE-102)  monitor  the  time  resolved  neutron  pulse  and  the  hard  X-ray  pulse 
emission  from  the  pinch  (X-ray  energy  >30  keV;  detector/monitoring  system  rise-time 
less  than  5  ns  ). 

We  note  that  the  non-optimized  APF-200  has  a  corresponding  distribution  over  a 
much  wider  interval  0. 1  <  Y/Y^^  <  14. 1,  with  Ygv( APF-200)  ~  Y,v(APF-50)/10. 

Fig.  6  reports  the  ion  spectrum  at  80°  from  a  Thomson  (parabola)  spectrometer. 
Differential  pumping  holds  a  10-^  Torr  vacuum  in  the  spectrometer  chamber,  which  is 
separated  from  the  discharge  chamber  by  a  150  -pm-diameter  pinhole  from  which  the  ion 
beam  is  extracted.  Fig.  7  reports  the  spectrometer  arrangement  with  respect  to  the  APF 
electrodes.  Other  characteristics  of  APF-50  were  extensively  described  in  previous 
reports  (  as  Stevens  Tech  Report  of  February  28,  1994,  submitted  to  the  AFOSR  under 
Contract  No.  AFOSR  -91-0158  and  previous  reports  under  the  same  Research  Contract). 

The  relatively-low  Y  values  for  Configuration  1  reported  in  Appendix  1 ,  as 
compared  to  the  corresponding  Y  values  for  Configuration  2,  at  a  lower  energy  W 
reported  in  this  Section,  indicate  the  difference  made  by  the  KE  component:  No  KE  was 
used  in  Configuration  1  (i.e.,  the  system  was  used  as  a  conventional  type  of  PF  machine, 
rather  than  an  APF  machine).  The  increase  by  a  factor  >  3  to  5  of  Y^^,  in  Configuration  2 
of  APF-50  is  a  direct  effect  of  the  introduction  of  KE  to  form  an  APF  from  a  PF. 

We  note  that  the  MAX  (Y)  value  is  less  affected  by  the  KE  than  Yg^  .  MAX(Y) 
has  essentially  the  same  value  for  both  APF  and  PF  mode  of  operation,  within  the 
uncertainty  (-20%)  of  the  silver  counter  determinations.  This  point  is  critical  to 
understand  the  role  of  the  interelectrode  (plasma)  current  density  and  its  distribution  in 
the  energy  densification  process  of  APF  and  PF  machines.  The  KE  reduces  substantially 
the  fluctuation  of  Y  from  pulse  to  pulse,  with  a  steady  shift  of  all  Y  values  toward 
MAX(Y). 
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Fig.  1;  APF-200  general  view. 
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Fig.  2.  Rogowski-coil  signal  dl(t)/dt  (upper  side  1)  and  current  signal  I(t)  V(t')  dt' 
from  a  self-integrating  capacitive  voltage  divider,  generated  in  the  same  two- 
module  APF  discharge  (25  kV,  Y  ~  YjJ.  The  capacitive  voltage  divider  in  this 
mode  of  operation,  provides  a  smoothed  signal,  independent  of  the  Rogowski-coil 
signal  for  measuring  one  and  the  the  same  electrical  parameter  I(t)  of  the  system. 


Fig.  3;  Distribution,from  a  run  of-  261  discharges  under  identical  conditions,  of  the 

number  of  APF-200  discharges  with  yield  between  and  YA^g^g  +  0.2  ,  at 

Vfl  =  25  kV  (W  =  35  kJ,  p  =  5.5  Torr  of  Dj).  Under  the  two-module  limitation  the 
APF  is  far  from  optimized  conditions  with  the  chosen  R,  r,  Lq,  1^^^.  ^o,  W,  but  it  it 
still  capable,  because  of  the  new  design  with  the  field  emitter  KE,  of  generating 
good  pinch  compression,  i.e. ,  to  have  in  typical  runs  of  100  to  300  discharges,  a 
Y,y  far  above  background,  by  a  factor  >  100  and  Max  (Y)  >  5  x  10’  2.5  Mev 
neutrons. 
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Fig.  4;  APF-50  general  view.  Configuration  1  at  right,  Configuration  2  at  left. 


Fig.  6:  Deuterium  ion  energy  spectrum  at  SO®.  Multiple  exposure  of  the  CR-39  target 
inside  the  Thomson  spectrometer,  on  which  the  ion  tracks  are  etched  in  a  NaOH 
solution  at  70  “C  after  exposure  (18  discharges  were  used  for  the  formation  of  the 
B-axis,  after  turning  off  the  ~  4  kV/cm  electric  field  between  the  compact 
spectrometer  magnetic  poles;  22  discharges  were  used  for  the  formation  of  the 
parabolic  ion  track  patterns  on  the  same  target. 
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MAX(Y)  is  a  figure  of  merit  more  closely  related  to  I  than  to  any  other  parameter 
of  the  optimized  machine  at  a  fixed  value  of  W,  and  corresponds  to  an  optimal  current 
density  distribution  between  the  electrodes.  This  optimal  current  distribution  is  more 
frequently  achieved  with  KE  than  without  KE. 

Other  tasks  of  the  program  directly  related  to  the  variations  of  the  physical 
characteristics  of  the  self-field-dominated  fragments  of  the  current  sheet,  were  supposed 
to  be  completed  in  the  originally-planned  second  and  third  year  of  the  research  program. 


Fig.  7  Thomson  spectrometer  arrangement  for  the  80°  measurements  in  APF-50. 
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Appendix  I :  W  >  20  kJ  PF  Machine  Specifications 

The  imploding  current  sheet  is  generated  by  a  plasma  focus  machine  (PF)  with  Mather's 
geometry  of  the  electrodes.  The  quality  of  the  pinch  and  of  the  radially-converging  plasma 
sheet,  with  a  maximum  current  density  between  the  electrodes,  is  conveniently  assessed 
from  magnetic  probe  data  and  also  via  the  neutron  yield  per  pulse  from  D  +  D  reactions. 
As  a  rule,  in  a  series  of  test  discharges  in  pure  deuterium  or  in  a  mixture  containing 
deuterium,  the  highest  neutron  yield  is  measured  when  the  interelectrode  current  has  a 
sharply-peaked  distribution  rather  than  a  diffused  space  distribution.  Neutron  yield 
measurements  for  a  pure-deuterium-filled  discharge  chamber  (Fig.  4(a)  of  Chapter  ^)  are 
reported  in  Table  I  from  each  of  two  series  of  >  50  discharges  at  a  relatively-Iow  voltage. 
Under  the  chosen  conditions,  the  occurrence  of  a  multiplicity  of  current  sheets  and  of 
related  complications  is  minimized.  Detailed  information  on  the  current  distribution  in  the 
sheet  and  the  sheet  velocity  v  are  obtained  by  magnetic  probe  measurements,  as  reported 
in-Fig.  1  of  this  Appendix.  The  first  of  two  probes,  (P-1  ),  is  always  located  just  above 
the  edge  of  the  hollow  center  electrode(anode),  i.e.,  at  a  radial  distance  from  the  electrode 
axis  equal  to  the  anode  outer  radius  ro=  26  mm.  The  second  probe,  (  P-2  ),  is  at  a  radial 
distance  of  1  cm  from  (P-1  ),  inside,  or  in  other  tests  outside,  the  anode  edge,  depending 
on-the  position  at  which  we  want  to  determine  v(r,  z).  Both  probes  are  on  the  same  plane 
orthogonal  to  the  electrode  axis,  at  a  distance  Zp  =  1  cm  above  the  edge  of  the  anode.  The 
outer  electrode  (cathode)  has  inner  radius  R=64  mm  (38  mm  interelectrode  gap).  The 
inductance  of  the  powering  capacitor  bank  of  total  capacitance  C  =  60  pp  (  two  modules, 
each  with  5  Aerovox  capacitors  of  30  nFI,  6  pF,  10.8  kJ  at  60  kV)  is  3  nH.  One  rail  gap 
switch  of  inductance  20  nFI  is  used  for  each  modules  (Maxwell  model  #  40200,  max. 
current  750  kA,  max.  voltage  100  kV;  jitter  <  2  ns).  The  two  switches  are  closed  by  a 
voltage  pulse  from  a  lOO-kV-multiple-trigger  generator  (Maxwell  mod.  #40151,  10 
outlets,  rise  time  <  1  ns,  jitter  <  10  ns),  activated  by  a  master  1 00-kV-trigger  (Maxwell 
mod.  #  40295  jitter  <  7  ns;  rise  time  <  0.2  ns). 


Table  1.  Leading  figures  of  merit  (peak  electrode  current  4^,  energy  input  W,  neutron 
yield  Y  per  pulse)  of  driver,  a  PF  with  max.  W -100  kJ  at  60  kV.  Period  9. 0+/-0. 1  ps. 


1  capacitor  bank  | 

neutron  yield  | 

neutron  yield 

1  D2  filling  1 

1  powering  energy  & 

per  pulse 

per  pulse 

pressure  for  | 

1  peak  electrode  current 

1 

maximum  value 

mean  value 

maximum  yield  | 

1  W  (kJ),  4,  (MA) 

! _ 

Max  Y(D+D) 

Y(D+D) 

1 

1 

P  1 

1 

i 

1  13  kJ,  0.8  MA(21  kV) 

0.5  x  lO'o 

1.2  X  109 

1 

5  Torn  1 

1  18  kJ,  1.0  MA  (25  kV) 

1 _  1 

nrM  •  t  . 

1.6  X  1010 

_ 1 

5.2  X  109 

6  Tort  1 

J  1 

The  inductance,  up  to  the  breech  end  of  the  electrodes,  including  power  transmission 
plates,  is  =  14  nFI.  The  inductance  of  the  interelectrode  gap,  from  breech  to  muzzle 
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side  of  the  anode  end  is  —  (iiy27t)  /„  ln(R/ro)  ^  (anode  length  4  =  145  mm)  The 
circuit  total  inductance  under  static  conditions  is  Lq  ~  +  L,  =  40  nH.  From  the 

measured  period  T  =  9.0'  +/-  0.1  |is  ~  27t  (  LC  of  the  complete  circuit  (powering 
system  and  electrodes,  with  the  moving  plasma  sheet  in  the  interelectrode  gap),  the 
corresponding  value  of  the  total  inductance  under  dynamic  conditions  is  L  =  34  +/-  1  nH. 
The  dynamic  inductance  accounts  for  the  resistivity  R  -  dL/dt  =  v  >  23  mOhm  to 
46  mOhm,  where  v  >4x10^  cm/s  to  10^  cm/s  is  the  measured  current-sheet  velocity  after 
rolling  off  the  interelectrode  gap,  till  the  beginning  of  the  radial  compression  phase 
( if  R  >  0,  then  L  -  [  1/2C  +  (1/40^  -  {2n/Ty  RV4  )i/2  ]  <  j. 

The  two  reversals  in  each  of  the  P-1  signals  in  Fig.  1(a)  and  (b)  indicate  that  three  current 
sheets  CS,,  CS2,  CS3  are  sweeping  the  magnetic  probe  in  these  discharges.  The  neutron 

yield  depends  on  specific  characteristics  of  the  layered  structure  of  the  current  distribution. 
Specifically,  Y  is  strongly  dependent  on  the  time  delay  At^  between  the  leading  CS,  and 
the  trailing  CS3  and,  to  a  lesser  extent,  on  the  delay  At,  between  CS,  and  CSj,  in 
agreement  with  the  equation  Y  r  Y  x[  a  (At,)-"  +  p  (At,)-^  ],  with  n  >  2,  and  a  «  P  . 
Y(W,p)  is  the  Y  mean  value  from  all  discharges  with  the  same  W,  p.  The  major  number  of 
discharges  has  a  sufficiently-high  value  of  At,  (as  those  of  Fig.  1),  so  that  the  first  a-term 
can  be  neglected.  By  expressing  At  in  nanosec  we  have  P(W,p)  =  1.6x10’  ns^.  Further 
insight  on  the  equation  relating  Y  and  the  At '  s  is  gained  from  its  derivation,  similar  to  that 
of  Hariline  and  Ratliff  equations  for  neural  networks  [I].  In  our  case,  the  i-th  response 
function  is  the  i-th  contribution  to  the  neutron  yield  from  the  i-th  imploding  current  sheet; 
the  input  intensity  (~  dl/dt)  is  determined  from  the  corresponding  probe  signal  of  the 
sheet  wliich  carries  a  fraction  I/I  of  the  total  electrode  current  I~  1  MA.  Additional 
informan'on  on  the  Y,  (At)s  correlation  coefficients  is  reported  in  Ref  [2]  of  thisAppendix. 


l.se  V/div  8.00  V  200  ns/div  -508.0  ns  1.50  V/dIv  0.00  V  208  ns/'dtv  -508.0  ns 


Fig.  1.  P-1  magnetic  probe  signals  frorn  tu'o  discharges  under  the  same  conditions  (W  =  18  kJ. 

p  =  6  Torr  of  pure  D2).  Y  =  0.35  Y  in  (a);  Y  =  0.9  Y  in  (b). 
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Abstract 

Our  experiments  have  established  strong  azimuthal  anisotropy  of  fusion  reaction  products  and  fast  ions  ejected  from  the 
plasma  during  focused  discharges  in  deuterium  or  a  mixture  of  deuterium  with  ^He  and/or  The  anisotropy  is  explained  in 
terms  of  localized  plasma  domains  with  enhanced  nuclear  reactivity  (emitting  isotropically),  displaced  from  the  axial  pinch 
region  by  up  to  5.5  cm.  The  displaced  domains  generate  and  trap  over  80%  of  the  observed  population  of  fast  ions  (E{>2  MeV), 
produce  the  bulk  of  the  '^0(d,  n)  reactions  {E^'^2  MeV)  in  the  plasma,  and  have  atomic  density  7x  10^*  cm“^.  Most  of 
the  D  +  D  and  ^He  +  D  fusion  yield,  which  requires  low  ion  energy  (£i<0.2  MeV),  comes  from  the  region  close  to  the  plasma- 
focus  symmetry  axis. 


1.  Introduction 

Plasma-focus  (PF)  discharges  are  well-known 
sources  of  extremely  high  neutron  and  high-energy- 
ion  yields.  Our  PF  device  has  a  Mather-type  elec¬ 
trode  geometry,  modified  with  a  field  distortion  ele¬ 
ment.  The  stored  energy  of  the  powering  capacitor 
bank  is  7  kJ  (at  V—  17  kV).  In  one  discharge,  the 
PF  machine  filled  with  deuterium  at  5  Torr  acceler¬ 
ates,  traps,  and  ejects  10'^  D"^  ions  (0.5<£'d<10 
MeV)  and  simultaneously  produces  10^  neutrons 
from  D(d,  n)  fusion  [1,2].  When  the  filling  gas  of 
the  PF  chamber  is  a  mixture  of  deuterium  (or  ^He) 
and  a  high-Z  element  ('^C,  ‘^O)  we  observe  a 

variety  of  (high-Z  nucleus)  +  (low-Z  nucleus)  nu¬ 
clear  reactions  with  apparent  reaction  thresholds  of 
0.5  MeV<£,f^<4  MeV  and  with  reaction  yields  of 


^  Work  was  supported  in  part  by  AFOSR,  contract  AFOSR-92- 
0158. 


I  O’- 10^  [2,3].  The  population  of  high-energy  ions 
can  be  linked  to  the  explosive  decay  of  the  magnetic 
structures  of  organized  plasma  filaments  [4,5].  The 
plasma  domains  of  diam.  <0,5  mm  with  high  neu¬ 
tron  and  fast-ion  emission  [  6-9  ]  have  been  observed 
at  a  distance  of  up  to  2  cm  from  the  pinch  region 
[8,10].  Also,  measurements  of  radioactivity  induced 
by  beams  on  end-on  targets  show  a  displacement 
from  the  z-axis  of  the  radioactive  spot  [11-13],  The 
displacement  of  the  plasma  domains  characterized  by 
high-intensity  ion  emission  will  induce  azimuthal 
emission  anisotropy  if  the  detectors  are  positioned  at 
distances  greater  than  but  comparable  to  the  domain 
displacement.  An  indication  of  the  appearance  of 
azimuthal  asymmetry  was  reported  in  Ref  [14], 
which  shows  that  the  deuteron-beam  intensity  ratio 
/(-4O°^0^O)°)//(O"^6^4O°)^2  for  £^>2.4 
MeV  and  approaches  unity  for£d<0.4  MeV. 

The  aim  of  this  work  is  to  prove  and  explain  the 
existence  of  the  pronounced  structure  in  the  azimu- 


0375-9601 /94/$07. 00  ©  1994  Elsevier  Science  B.V.  All  rights  reserved 
5'5©/  0375-960I  (  94)00533-U 


18 


J.S.  Brzosko  et  ai  /  Physics  Letters  A  192  (1994)  250-257 


251 


lhal  distribution  of  the  ion  and  fusion-product 
emission. 


2.  Experimental 

The  space  distribution  of  plasma  sources  emitting 
fast  ions  and  reaction  products  is  monitored  (see  Fig. 
1 )  with  fifteen  thin-wall  Geiger-Mulier  counters  ( 1 5- 
GM).  The  counters  are  parallel  to  the  PF  symmetry 
axis  in  a  12.6  cm  diameter  circle.  The  chamber  walls 
and  the  stopping  disc  at  the  front  end  of  the  chamber 
are  clad  with  tantalum,  which  is  thick  enough  to  stop 
ions  with  energy  ^  20  MeV.  The  counters  measure 
the  radioactivities  from  reaction  products  deposited 
on,  or  produced  in,  the  metal  foils  wrapping  the 
counters.  Each  wrapped  detector  covers  /i6=  1 T  and 
25°<(;^j<60'’;  where  6  is  the  azimuthal  angle  and  tp  is 
the  angle  between  the  PF  symmetry  axis  and  the  po¬ 
sition  vector  originating  from  the  center  of  the  anode 
front  end.  Radioactivities  measured  simultaneously 
for  all  counters  are  used  to  reconstruct  the  spatial  po¬ 
sition  (s)  (R,  6)  and  intrinsic  radioactivity  (ies)  Aq 
of  the  emitting  source (s)  (see  Fig.  Ic). 

Within  60  s  of  each  discharge,  the  PF  chamber  is 
evacuated  and  the  15-GM  system  begins  recording 
the  radioactivity  decay  curves,  separately  for  each  GM 
counter.  The  decay  curves  are  stored  by  the  computer 
and  the  initial  activities  {A),  defined  at  pinch  instant 
and  corresponding  half-lives  (rj/2)  of  the  decaying 
nuclei,  are  determined  by  standard  numerical  meth¬ 
ods  ( least-squares  method  with  /2  and  A  as  free  pa¬ 
rameters).  From  the  determined  r,/2,  we  unambig¬ 
uously  define  the  emitting  radioisotope.  The  final 
reaction  product,  ion  energy  spectrum  [  1  ],  and 
isotope  composition  in  the  PF  chamber  define  the 
breeding  nuclear  reaction.  The  rj/2  values  obtained 
from  this  procedure  agree  with  the  commonly  ac¬ 
cepted  values  [15]  to  within  a  few  percent. 

During  the  experiment,  the  detector  efficiencies  are 
stable  to  within  1%  and  the  cross-talk  between 
counters  (defined  experimentally)  does  not  change 
significantly  the  azimuthal  distribution  of/f. 

The  neutron  yield  per  discharge  from  D(d, 
n)^He  fusion  reactions  in  the  plasma  was  measured 
with  two  identical  Ag-activalion  counters  (Ag-GM) 
of  LANL  design  [16].  The  Ag-GM  counters  were  po¬ 
sitioned  symmetrically  to  the  PF  electrode  axis  to 


Fig  1.  Experimental  set-up.  (a)  Axial  cross-section  of  the  plasma- 
focus  chamber  (1)  PF  coaxial  electrodes;  (2)  insulator,  (3) 
chamber  walls  clad  with  tantalum;  (4 )  one  of  the  1 5-GM  counters 
(wrapped  with  thin  AJ  or  Cu)  set  parallel  to  the  PF  symmetry 
axis;  (5)  Ta-stopping  disc,  (b)  Section  in  the  plane  orthogonal 
to  the  PF-axis.  GM-detectors  are  identified  by  an  azimuthal  an¬ 
gle  {6)  assigned  to  the  detector  anode-wire.  The  star  marks  one 
displaced  PDER.  (c)  Expected  count  distribution  (see  Eqs.  (2) 
and  (3))  for  a  PDER  disintegrating  in  the  position  shown  in  (b). 

measure  the  azimuthal  anisotropy  of  the  neutron 
emission. 

We  also  routinely  recorded  the  time  variation  of 
dl/dt  (/is  the  electrode  current  measured  by  the  Ro- 
gowski  coil)  and  of  dY^/dt  (detected  by  an  NE102 
plastic  scintillator)  using  two  channels  of  a  digitized 
scope.  The  time  resolution  of  these  measurements  is 
about  5  ns.  These  data  confirm  the  well-known  fact 
[17-19]  that  the  high  nuclear  reactivity  phase  of  the 
plasma  occurs  after  the  maximum  axial  compression 
of  the  pinch.  In  our  device  the  delay  between  the  peak 
values  of  \dl/dt\  and  that  of  |d7„/d/|  is  typically 
zl/=  100  ns.  We  also  found  that  dTn/d/  and  dl/dt  data 
are  not  clearly  related  to  the  observed  azimuthal 
asymmetries  of  the  ion  and  fusion-product  emission. 
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3.  Trends  in  experimental  data 

The  experimental  program  is  based  on  three  com- 
plementar)’  experiments,  each  with  a  different  gas 
mixture  filling  the  discharge  chamber  and  different 
material  wrapping  the  GM-counters.  For  a  complete 
specification  of  the  three  experiments  and  of  the  en¬ 
ergy  threshold  of  the  observed  nuclear  reactions,  see 
Table  I. 

Measurement  of  the  azimuthal  distribution  of  the 
reaction  products  in  the  plasma  and  the  products  of 
reactions  induced  by  the  D*^  beam  in  the  external 
targets  (wrapping  foils)  indicates  the  following: 

(a)  Only  one  half-life  decay  curve  is  observed  in  ex¬ 
periment  1  (D2  filling,  Al  wrapping).  The  corre¬ 
sponding  activity  is  identified  as  the  product  of  the 
reactions  induced  in  the  aluminium  wrapping  foil  by 
the  D'^  beam,  i.e.,  ^’Al (d,  p)^®AL 

(b)  A  double  half-life  decay  curve  is  observed  in  ex¬ 
periments  2  and  3.  In  experiment  2,  one  half-life  re¬ 
sults  from  the  decay  of  products  of  the  *®0(d,  n)'’^F 
reaction  in  the  plasma;  the  second,  from  the  decay  of 
products  of  the  ^^Cu(d,  p)^^Cu  reactions  in  the 
wrapping  foil.  In  experiment  3,  one  half-life  results 


from  the  decay  of  ^^Zn  nuclei  produced  in  the  chain 
of  two  reactions — the  ^He(d,  p)'*He  reaction 
(£p=  14.7  MeV)  in  the  plasma  followed  by  ®^Cu(p, 
n)^^Zn  in  the  wrapping  foil;  the  second,  from  the  de¬ 
cay  of  products  of  ^^Cu(d,  p)^^Cu  reactions  in  the 
wrapping  foil  (as  in  experiment  2).  Consequently,  we 
simultaneously  defined  two  azimuthal  angular  distri¬ 
butions — one  for  fusion  products,  the  other  for  fast 
deuterons. 

(c)  The  observed  distributions  have  a  strong  azi¬ 
muthal  structure  in  more  than  80%  of  the  discharges, 
with  up  to  three  distinct  maxima.  Examples  of  the 
observed  distributions  are  shown  in  Figs.  2,  3,  and  4. 

(d)  There  is  a  strong  correlation  between  the  azi¬ 
muthal  distributions  measured  in  the  same  discharge 
for  radioactivities  induced  by  fusion  processes  and 
by  D*^  beams  on  the  targets  (wrapping  foils);  see  Figs. 
3  and  4. 

(e)  The  positions  of  the  maxima  of  the  detected  ra¬ 
dioactivities  are  randomly  distributed  along  the  azi¬ 
muthal  coordinate  6. 

(f)  The  azimuthal  neutron  anisotropy  is  ran¬ 
domly  distributed  between  —0.1  and  +0.1,  and  is 
defined  as 


Table  1 

Specification  of  the  experiments 


Gas 

composition  * 

Foil 

[mm] 

Monitoring  reaction 

beam 

plasma  fusion 

reaction 

[MeV] 

reaction 

[MeV] 

expenmenl  I 

D2(100%) 

Al 

^’AI(d,p)^'Al 

1.8 

D(d.n)5He'’ 

0.005 

0.254 

7,;2  =  2.4  m 

experiment  2 

D2(90%) +02(10%) 

Cu 

‘^Cu(d,p)“Cu 

2.0 

D(d,n)^He'’ 

0.005 

0.127 

7'i/2  =  5.1  m 

'‘0(d.n)'’F' 

2.0 

Ti,i=\.\  m 

experiment  3 

0.(50%) +^He(50%) 

Cu 

‘-Cu(d,p)'’'’Cu 

2.0 

D(d,n)’He'> 

0.005 

0.127 

7'i/2  =  5.I  m 

^He(d,p)^He  then" 

0.1 

“Cu(p,n)"’Zn 

5.0 

r,/2  =  38.1  m 

’  The  chamber  is  filled  with  3.3  x  1  O'"  atoms/cm''  (equivalent  to  5  TorrorD2).  Numbers  in  brackets  refer  to  atomic  abundance. 
Azimuthal  anisotropy  {A^)  of  neutron  emission  was  measured  only. 

'"F  ions  (£,  ca.  0.6  MeV^ )  are  implanted  in  Cu-foil  wrapping  GM-counlcrs. 

14,7  MeV  proions  from  ^He(d,p)^He  reaction  in  the  plasma  were  delected  via  the  (p,n)®'^Zn  reaction  induced  in  Cu-foil  wrapping 
GM  counters. 

is  the  apparent  reaction  threshold  defined  as  the  ion  energy  above  which  90%  of  reactions  occur  (for  ion  spectrum  see  Ref 

i  1  ]  )• 
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(a)  ( 

■ 

y,  ^ 

c  0  250 


Fig.  2.  Examples  of  azimuthal  distributions  of  radioactivity 
produced  in  the  ^’Al(d,  p)^®Al  reaction  by  deuterons  with  energy 
£’d>  2  MeV  (experiment  1 ).  Full  squares  with  error  bars  (stan¬ 
dard  deviation  of  measured  radioactivities)  are  the  experimental 
points.  Histograms  are  the  best  fit  of^i,  0^  and  R,  to  the  function 
representing  anisotropy  as  the  result  of  different  azimuthal  and 
radial  positions  of  emitting  PDERs  (linear  combination  of  terms 
as  in  Eq.  (2) ).  (a)  An  example  of  a  PDER  on  impact  with  an  Al 
foil  wrapping  a  GM-counter  at  0=264®  (/I  =  61%);  remaining 
yield  is  divided  into  a  PDER  at  .^2  =  5.0  cm  (02=  16®,  ^^12=  16%) 
and  a  PDER  on  the  PF-axis  (^  =  23%).  (b)  A  single  PDER  at 
R  =  4.S  cm,  d-  -48“  carrying  100%  of  the  D"^  beam,  (c)  Two 
displaced  PDER  sources:  at  /?|  =5. 1  cm  (0,  =  13®, /I,  =  52%)  and 
/?2  =  4.8  cm  (02  =  86,  /42  =  28%)  and  an  axial  PDER  (/43  =  20%). 
(d)  An  example  of  multiple  PDER  emission:  /?,=4.2  cm 
(0,=  13O®,  /I,  =  24%),  y?2  =  4.2  cm  (02  =  230®,  A2  =  24%)  and 
/?3  =  5.2  cm  (03=333®,  ^13  =  32%)  R^==5.2  cm  (0^  =  381“, 
/l4  =  20%). 


o 

< 

't) 


Azimuthal  Angle  -  0° 

Fig.  3.  An  example  of  comparison  between  fusion-  and  beam- 
generated  radioactivities  as  recorded  from  the  same  PF  discharge 
by  the  1 5-GM  system.  The  PF  machine  is  filled  with  D2  (2.5  Torr) 
and  ^He  (5  Torr).  Experiment  2:  (a)  PDERs  (beam)  observed 
at  R^  =  5,5  cm  (0,  =  6O®.  .4,  =  12%)  and  /?2  =  4.0  cm  (02  =  241®, 
.12  =  74%)  and  a  central  PDER  of/l3=  14%;  (b)  a  PDER  (^He  +  D 
fusion )  measured  by  radioactivity  induced  by  14.7  MeV  protons 
and  observed  at  /?,  =  5.5  cm  (0,  =40®,  /I,  =  6%)  and  /?2  =  4.0  cm 
( 02  =  24 1 ®,  A2  =  32% )  and  an  axial  PDER  of /I3  =  62%. 


/ln(/=23  cm)  =  1  -  r,(0=9O'’)/rj(9=:27O°)  . 

(la) 

The  neutron  emission  anisotropy  is  measured  by 


Fig.  4.  Similar  to  Fig.  3  but  the  PF  machine  is  filled  with  D2  (4.5 
Torr)  and  O2  (0,5  Torr).  Experiment  3:  (a)  PDERs  (beam)  ob¬ 
served  at  /?,=4.6  cm  (0,  =34®,  ,^,=69%)  and  /?2  =  3.3  cm 
(02=235®, /l2  =  3l%);  (b)  a  PDER  C^O  +  D  fusion)  measured 
by  ‘^F  radioactive  nuclei  implanted  into  the  Cu  wrapping  foil 
and  observed  at i?,  =  4.5  cm  (0,  =  33®,v4,  =  95%)  and  7^2  =  3. 3  cm 
(02=  192®,/l2  =  3%). 


the  two  Ag-GM  neutron  counters,  which  are  posi¬ 
tioned  1=23  cm  from  the  PF  axis,  on  opposite  sides 
of  it. 

(g)  The  linear  correlation  coefficient,  p  between 
neutron  anisotropy  ^n(^=  23  cm)  and  (or  ^^Cu) 
radioactivity  anisotropy  (.4^)  for  21  discharges  is 
p{A^\  A^)  =  0A  with  a  confidence  level  of  85%  (con¬ 
fidence  level  defined  as  in  Ref  [20]).  A^.  is  esti¬ 
mated  from  the  activity  measured  by  GM  counters 
covering  the  same  azimuthal  angles  as  covered  by  the 
neutron  counters  and  is  defined  as  follows: 

A^  =  \~A{4S^<e<\20^)/A{2W<e<3\2°)  . 

(lb) 

Only  those  discharges  where  the  sum  of  the  ra¬ 
dioactivities  measured  by  the  counters  used  to  define 
Ad  is  at  least  40%  of  the  total  radioactivity  measured 
by  all  the  15-GM  counters  are  taken  for  analysis.  A 
correlation  test  {An,  A^ )  with  the  neutron  counters  at 
/=98cm  does  not  reveal  any  significant  correlation. 

Measurements  were  made  of  the  dependence  of 
neutron  asymmetry  An  on  the  Ag-GM  counter  dis¬ 
tance.  Two  anisotropy  tests  were  conducted  at  1=23 
cm  (100  discharges)  and  at  /=98  cm  (50  dis¬ 
charges).  We  also  compared  the  detector  readings  at 
/=40  cm  (50  discharges)  and  at  azimuthal  positions 
as  close  as  permitted  by  the  experimental  conditions 
(within  A0=20°)  to  check  the  detector  stability.  For 
/^n(/=98)  and  An{A6=20°),  the  results  are  well  rep¬ 
resented  by  a  normal  distribution  (with  standard  de¬ 
viation  (7),  which  is  governed  by  statistical  errors  of 
the  counter  readings,  e.g.,  a/AAn=  1.  On  the  other 
hand,  for  /!„(/=  23  cm),  30%  of  the  events  differs 
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from  the  average  value  by  more  than  one  standard 
deviation  ( cj)  and  a/AA^  >  3.  The  latter  is  due  to  the 
radial  and  azimuthal  displacement  of  the  neutron 
source.  Such  a  displacement  afTects  the  readings  of 
the  detectors  at  small  values  of  distance  /  and  can  be 
neglected  for  detectors  far  from  the  z-axis. 


4.  Data  interpretation 

The  analysis  of  the  azimuthal  anisotropy  observed 
in  the  experiments  assumes  that  any  part  of  the 
plasma  volume  (defined  by  the  1 5-GM  counters,  the 
tantalum  stopper  in  front,  and  the  plane  determined 
b)’  the  front-end  of  the  external  electrode)  may  iso- 
tropicalh  emit  ions  and  fusion  reaction  products. 
Consequently,  the  anisotropy  will  be  the  result  of  dis¬ 
placed  (from  the  PF-axis),  localized  plasma  do¬ 
mains  of  enhanced  nuclear  reactivity  (referred  to 
henceforth  as  PDER(s) ). 

The  alternative  assumption  of  an  anisotropic  (azi¬ 
muthal)  emission  of  ions  and  fusion  products 
from  the  central  pinch  column  is  rejected  because  it 
IS  unable  to  consistently  explain  all  the  experimental 
data  in  this  work.  We  refer  specifically  to  the  data 
quoted  in  points  (d)  and  (g)  above. 

Our  experimental  system  is  sensitive  to  the  radial 
positionfs)  of  emitting  plasma  domain(s)  {R)  and 
its  azimuthal  coordinate(s)  (Oq).  but  the  experi¬ 
ments  do  not  allow  us  to  define  the  axial  position (s) 
of  the  domain  (s).  As  a  consequence,  a  PDER  emit¬ 
ting  during  impact  on  the  Ta-stopper  cannot  be  dis¬ 
tinguished  from  a  PDER  in  radial  flight  through  the 
plasma  volume.  To  construct  a  map  of  PDERs  that 
would  explain  the  observ^ed  anisotropy  one  can  use 
the  least-squares  method  [20]  to  compare  radioac¬ 
tivities  imposed  by  the  best  fit  of  source  distributions 
(D^  beams  and/or  fusion  products)  in  the  plasma 
with  the  measured  radioactivities.  The  tests  were  done 
for  two  space  representations  of  radiating  plasmas. 

In  the  first  representation,  the  cylinder  of  plasma 
normal  to  the  PF  axes  ( I  cm  thick,  with  base  at  z>  5 
^  in )  is  divided  into  thirteen  areas  and  the  thirteen 
intrinsic  radioactivities  (A^^)  can  be  fitted.  The  azir 
muthal  division  of  the  plasma  source  is  represented 
by  a  polar  surface  containing  a  central  disc  plus  six 
radial  and  two  azimuthal  grids  such  that  all  thirteen 
iiave  equal  surface  area.  The  second  model  considers 


point-like  plasma  sources  (axial  and  up  to  four  ra¬ 
dial)  with  Aq,  R  and  coordinates  as  parameters  to 
be  fitted.  As  a  special  case  of  the  second  model  one 
may  consider  that  the  emission  from  the  plasma  can 
occur  with  plasma  lumps  impacting  on  the  surfaces 
of  the  foils  wrapping  the  1  5-GM  counters.  In  any  case, 
the  optimization  procedure  (least-squares  method) 
fits  an  analytical  function  A^y^{0)  to  the  experimen¬ 
tally  determined  radioactivities >4„p(^).  Any  plasma 
emitting  volume  is  represented  by  superposition  of 
emitting  point-iike  sources. 

For  one  point-like  PDER  at  the  (R,  do)  location, 
the  distribution  of  radioactivity is  written  as 

A,^(e)=AZQ(R,eo,e)  (2) 

and 

2 

Q{R,  6q,  6)--  -arctg(<3/?r“‘(<3^4-/5^  +  r‘)  ^ 

(3a) 

r"  =  i?--f-L"~2RL  cos(/9o  “^)  ■  (3b) 

The  emitting  point-like  plasma  domain  is  localized 
at  a  distance  R  from  the  PF-axis  at  z=a/2  and  at  an 
azimuthal  angle  6^:  L  — 6.3  cm  is  the  distance  from 
the  GM  axis  to  the  z-axis;  ^=10  cm  and  /7  =  0.8  cm 
are  the  active  length  and  radius  of  the  GM  counters. 
Fig.  Ic  shows  the  A^^{6)  distribution  that  would  be 
observ^ed  if  a  single,  point-like  PDER  is  in  the  posi¬ 
tion  marked  in  Fig.  lb. 

The/^  test  w’as  used  to  choose  between  a  represen¬ 
tation  of  plasma  sources  by  13  large  elements  (with 
constant  emissivity  within  each  area)  and  point-like 
sources  with  (R.  6)  coordinates  as  free  parameters. 
For  a  few  typical  experimental  data  sets,  the  best  fits 
were  done  and  yr/xh  was  defined  (see  Fig.  5).  yl^ 
corresponds  to  a  fit  with  a  confidence  level  of  50% 
and  is  tabulated  elsewhere  [20].  This  analysis  leads 
to  the  following  conclusions: 

(i)  all  the  considered  representations  of  plasma 
sources  give  a  similar  azimuthal  position  of  the  main 
emitting  plasma  source  but  a  different  angular  width 
of  the  source(s ); 

(ii)  the  highest  confidence  level  is  achieved  for  a 
small  number  of  point-iike  sources,  where  the  num¬ 
ber  of  observed  maxima  in  the  distribution  ( if  clearl  v 
separated)  is  equal  to  the  number  of  PDERs  dis¬ 
placed  more  than  r>  3  cm; 
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Fig.  5-  of  the  different  localized  plasma  domains,  x^lxl.^ 
is  plotted  as  a  function  of  the  number  of  emitting  domains  and 
inserts  show  the  best  fit,  e.g.,  for  the  case  when  x^  1x1.5  a  mini¬ 
mum.  The  following  plasma  representations  are  used:  the  point 
sources  ( □ );  the  lump  sources  on  impact  with  a  set  of  neighbor¬ 
ing  GM-counters  (A);  and  the  1 3  equal-volume  sources  {0)^X^ 
is  from  comparison  ofA^  and  A„p\  xh  is  the  theoretical  value  of 

that  would  be  obtained  for  a  fit  with  a  confidence  level  of  50%. 

(iii)  will  not  change  significantly  if  one  substi¬ 
tutes  point  sources  with  spheres  (centered  at  the  lo¬ 
cation  of  the  point  source)  having  radius  k  smaller 
than  the  distance  from  the  sphere  center  to  the  walls 
of  the  closest  GM  counter;  in  particular,  for 
2k<L—R,  x^  will  change  by  less  than  20%. 

For  further  discussions  we  will  use  the  point-like 
representation  of  PDERs  expressed  by  Eqs.  (2)  and 
(3). 

The  analysis  of  50  azimuthal  distributions  from 
experiments  1  and  2  shows  that  in  more  than  80%  of 
the  discharges  we  can  identify  at  least  one  (typically 
two)  high-intensity  PDERs  displaced  by  more  than 
3  cm  from  the  z-axis.  Out  of  all  the  displaced  PDERs, 
35%  disintegrates  on  impact  on  the  GM  wrapping  foil 
{R  =  55  mm±l  mm).  One  can  also  conclude  that 
PDERs  disintegrate  rather  far  from  the  PF-axis,  while 
PDER  ion  emission  yield  remains  quasi  independent 
(within  ±  50%)  on  its  disintegration  position.  It  was 
found  that  displaced  PDERs  systematically  emit  more 
than  80%  of  the  total  amount  of  the  ejected  fast-ion 
(Ei>2  MeV)  and  *^0(d,  n)‘’F  reaction  products 
(£th>2  MeV).  The  fits  indicate  that  in  10%  of  the 
discharges,  the  fast  ions  come  from  one  displaced 
PDER  only. 

The  data  of  experiment  3  (D2  +  ^He  filling)  reveal 
a  correlation  between  the  azimuthal  distributions  of 
the  radioactivity  from  ^^Cu  (monitoring  £^>2  MeV) 
and  from  ^^Zn  (monitoring  the  14  MeV  protons  from 
the  ^He(d,  p)'*He  fusion  induced  by  deuterons  with 
Ed  <  0.2  MeV ) ;  see  Fig.  3.  Numerical  analysis  of  these 


data  shows  the  same  PDER  positions  (R,  6q)  and 
substantial  differences  in  the  relative  yields  of  dis¬ 
placed  and  axial  sources  for  both  ^^Cu  and  ^^Zn  ra¬ 
dioisotopes  produced  by  high-energy  and  low-energy 
deuterons,  respectively.  Most  (70%-90%)  of  the  ^^Zn 
radioactivity  is  associated  with  an  axial  (symmetric ) 
emission  of  14  MeV  protons  and  only  10%-30%  is 
related  to  emission  from  displaced  PDERs,  while  the 
^^Cu  radioactivity  (from  the  same  discharge)  is  as¬ 
sociated  mainly  with  £^>2  MeV  deuterons  emitted 
by  displaced  PDERs.  Both  cases  (14  MeV  protons 
and  Ed>2  MeV  deuterons)  reveal  the  same  posi¬ 
tions  (Rj,  doi)  for  the  displaced  PDERs,  On  the  other 
hand,  the  axial  position  of  the  main  source  of  the 
^He+D  fusion  is  in  agreement  with  the  well  estab¬ 
lished  fact  of  the  axial  location  (pinch )  of  the  leading 
D  +  D  fusion  source  [21-24]. 

The  neutron  yield  distribution  between  the  axial 
plasma  column  and  displaced  PDERs  was  estimated 
from  a  comparison  of  the  measured  neutron  anisot¬ 
ropy  y4„(/=23  cm)  and  the  neutron  anisotropy 
/t„(PDER)  inferred  from  the  localization  {R^,  ^oi)  of 
PDERs  (found  from  the  distribution  of  radioiso¬ 
topes  induced  by  deuterons  of  Ed >2  MeV).  In  ad¬ 
dition  to  the  positions  and  relative  intensities  of  the 
PDERs,  we  assumed  that  in  each  discharge  80%  of 
the  neutron  emission  is  from  the  axial  pinch.  This  as¬ 
sumption  allows  a  good  reproduction  of  the  experi¬ 
mental  data  (see  Fig.  6). 

A  lower  bound  of  a  PDER  atomic  density  n  (time 


Fig.  6.  Comparison  of  the  measured  neutron  emission  anisotropy 
A^  (two  Ag-GM  placed  at  0,  =90’  and  at  02  =  270’  with  /=23  cm) 
with  the  estimated  y4o(PDER)  from  positions  and  intensities  of 
PDERs  (estimated  from  data  of  experiment  1 )  with  the  addi¬ 
tional  assumption  that  80%  of  the  neutrons  are  emitted  from  the 
PF-axis.  The  error  bars  have  their  origin  in  the  statistical  errors 
of  the  neutron  detector  counts. 
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averaged  during  the  nuclear  reactiviiy  phase)  can  be 
estimated  from  the  experimental  data  of  Fig.  4  by  us¬ 
ing  similar  approaches  to  those  used  in  Refs.  [2,3], 
For  this  estimate  we  have  to  compare  the  ratio  of  the 
total  radioactivity  induced  by  one  PDER  via  ’^0(d. 
n)'^F  reactions  in  the  plasma  with  the  total  radioac¬ 
tivity  induced  by  ejected  from  the  same  PDER. 
The  reaction  cross  sections,  the  spectra  of  trapped  and 
reacting  ions  [1],  and  the  PDER-Iifetime 

fraction  during  which  ion  acceleration  and  trapping 
occur  (  ~  5  ns  [  8  ] )  are  the  known  quantities  entering 
this  estimate.  From  the  data  shown  in  Fig.  4  w'e  found 
that  n  ^  7  X  1 0^  ‘  /cm  "  This  density  is  typical  of  all 
the  data  from  experiment  3.  It  also  agrees  with  the 
estimation  in  Ref.  [  2  ]  that  is  based  on  the  data  from 
space-integrated  fast-ion  fusion  and  fast-ion  emis¬ 
sion.  Such  a  large  plasma  density  in  the  plasma  focus 
discharge  (pinch  density  is  10 ’^/cm^)  was  deter¬ 
mined  in  the  past  via  CO2  laser  scattering  on  50-/im 
plasma  domains,  n>  lO^^/cm^  [25]  and  from  X-ray 
line  intensities  from  'Tot-spots”,  n>  10^‘/cm^  [26]. 

5.  Summan’ 

All  the  data  presented  in  this  paper  can  be  consis¬ 
tently  explained  by  displaced  plasma  domains  of  very 
high  density.  It  is  conceivable  that  the  displaced 
PDERs  and  hot-spots  (diam.<0.5  mm)  previously 
observed  in  the  vicinity  of  the  r-axis  by  pin-hole  im¬ 
aging  (D"^.  Ea>  1  MeV,  recorded  with  CR-39  track 
detectors)  [8,27]  have  the  same  origin  and  physical 
characteristics. 

The  following  conclusions  are  drawn: 

( 1  )  High-energy  ions  (Ed >  2  MeV )  and  the  prod¬ 
ucts  of  the  '^0(d,  n)^'^F  nuclear  reactions  (reaction 
threshold  E,h>  2  MeV)  in  the  plasma  are  both  emit¬ 
ted  from  the  same  PDER.  The  major  fraction 
(-'80%)  of  the  PDERs  is  displaced  from  the  PF 
symmetry’  axis  at  a  distance  greater  than  3  cm,  and 
about  35%  of  the  population  is  displaced  up  to  5.5 
cm  (a  limit  set  by  the  experimental  conditions). 
Often,  more  than  80%  of  the  deuterons  of  the  Ed>  2 
MeV  and  80%  of  the  reaction  yields  induced  by  these 
deuterons  in  the  plasma  are  generated  by  one  PDER 
displaced  far  from  the  z-axis.  The  remaining  20%  of 
the  emission  is  emitted  from  the  PF  symmetry  axis. 

( 2 )  The  displaced  PDERS  are  also  sources  of  D  -f  D 


and  D  +  ^He  fusion  reactions  induced  by  ions  of  en¬ 
ergy  Ed  ^0.2  MeV.  However,  only  10%-30%  of  the 
D  +  D  and  D  +  ^He  fusion  reactions  are  associated 
with  displaced  PDERs.  Most  of  the  D  +  D  and 
D  +  ^He  fusion  reactions  have  their  origins  in  the  near 
PF-axis  region,  as  reported  elsewhere  [21-24]. 

This  observation  introduces  spatial  differentiation 
between  plasma  domains  poor  in  high-energy  ions  and 
those  rich  in  high-energy  ions.  A  comparison  of  the 
present  results  with  those  of  Ref.  [2]  suggests  that 
displaced  and  axial  plasma  domains  have  different 
ion  energy  characteristics  but  quadratic  relations  of 
their  nuclear  reactivities. 

The  PDER  displacement  may  be  enhanced  by  the 
grounding  of  the  discharge  chamber  walls,  GM- 
counter  wrapping  foils,  and  the  front  stopper.  If  this 
is  so,  it  would  be  an  important  process  for  future  ap¬ 
plications  of  PF  machines. 
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